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Context. Light curves of the long-period Algols are known for their complex shape (asymmetry in the eclipse, light variations 
outside eclipse, changes from cycle-to-cycle), and their interpretation is not possible in the standard model of binary stars. 
Aims. Complex structures present in these active Algol systems could be studied with the eclipse-mapping method which was 
successfully applied to the new 7-color photometric observations in the Geneva system of W Cru, belonging to the isolated 
group of these active Algols. 

Methods. Several cycles of this long-period (198.5 days) eclipsing binary have been covered by observations. We have used a 
modified Rutten's approach to the eclipse-mapping. The optimization of the system's parameters and the recovery of the disk 
intensity distrubution are performed using a genetic algorithm (GA). 

Results. It is found that a primary (hot) component is hidden in the thick accretion disk which confirms previous findings. The 
Q mass of the primary component. Mi — 8.2 M0 indicates that it is a mid-B type star. The mass-losing component is filling 
^ ^ its critical lobe which, for the system's parameters, means it is a G-type supergiant with a mass M2 = 1.6 Mq. The disk is 
CZ3 I geometrically very extended and its outer radius is about 80% of the primary's critical lobe. A reconstructed image reveals 
a rather clumpy and nonuniform brightness distribution of an accretion disk rim in this almost edge-on seen system. This 
dumpiness accounts for light curve distortions and asymmetries, as well as for secular changes. 
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1. Introduction 

W Cru (HD 105998) is a member of the rather sparse 
group of strongly interacting binaries which are believed to 
be either in the first and rapid phase of the mass exchange 
between components, or have just recently gone it. While 
mass reversal has already happened there are still many 
signatures of the almost cataclysmic events that were tak- 
ing place. Studying the UV spectra of some long-period 
binaries taken with the spectrograph on-board lUE satel- 
lite Plavec & Koch (1978) have found numerous emission 
lines of highly-ionized species. It was immediately obvious 
that the underlying stellar photospheres of these binaries 
are not hot enough to excite such features, and therefore 
the accretion processes were invoked as an explanation. 
This sparse group of the long-period interacting binaries 
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comprises SX Cas, RX Cas, W Cru, W Ser, V367 Cyg, and 
(3 Lyr, and a few others less studied systems. All these bi- 
naries are well known for their active nature with large 
period changes, almost permanent Balmer line emission 
features, peculiar light curves, etc. (Plavec 1980). 

In the long-period Algol systems with P > 6 days, 
the primary component is rather small relative to the bi- 
nary separation and its Roche lobe, and mass transfer 
results in a classical accretion disk. Typically, this disk 
is a permanent structure characterized by strong double- 
peaked emission lines in the spectrum. Since the disk hides 
the mass-gaining component, which in these post mass- 
exchange binaries is more massive and hotter than its 
companion, we are encountering a paradox that the com- 
ponent expected to be more luminous is invisible. This 
effect considerably affect the light curves. 

W Cru is not an unique binary system which har- 
bors such a large accretion disk. Several other long-period 
semi-detached binary systems have been found for which 
light curve solution has been found only after an opti- 
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Fig. 1. The light and colour curves of W Cru in the 

Geneva photometric system. The light curves in G and 
Bl have been vertically shotted by 0.3 mag for clarity. 



cally and geometrically thick disk has been taken into 
account. Good examples are, RZ Oph (Zola (1991), and 
UU Cnc (Zola et al. 1994), with the periods of 262 and 
96 days, respectively. In all these long-period Algols the 
separation between the components is large, and hence, 
mass-losing components are supergiants, and disks which 
surround mass-accreting (hot and more massive) compo- 
nent typically extend several hundred solar radii in radial 
direction. 

It is well known that light curves of the 'active al- 
gols' cannot be solved with standard models that deal 
only with two (stellar) components. In this respect the 
puzzling case of /? Lyrae is well-known (c.f. Linnell 2000, 
Harmanec 2002, and references therein) but the situation 
with other binaries of this class is quite similar. A real 
break-through came with Huang's (1963) disk-model for 
the p Lyr, and its first application to V356 Sgr by Wilson 
& Caldwell (1978). Following this development Pavlovski 
& Kfiz (1985) made a computer code for the synthesis of 
light curves of close binary stars in which the influence 



of an accretion disk is taken into account. They applied 
it in an analysis of the light curves of SX Cas. Derived 
disk properties were in general accordance with multi- 
wavelength studies by Plavec et al. (1982). Therefore, 
subsequent studies have been performed on more difficult 
cases: RX Cas (Andersen, Pavlovski & Piirola 1989) and 
V367 Cyg (Pavlovski, Schneider & Akan 1992). 

In the meantime, some other attempts along the same 
disk-model but with more or less physics details and var- 
ious optimization techniques have been undertaken. In 
particular, Zola (1991) has implemented a disk-model in 
the widely used Wilson-Devinney program, while Daems 
(1998) introduced an analytical disk model. However, 
these studies do not take into account the fact that the 
accretion disk intensity distribution is generally asymmet- 
ric. This means that a more realistic model(s) is needed 
which is capable to account for erratic changes in the light 
curves, asymmetry in ascending and descending branches, 
unequal brightness of maxima, etc. In other words, a dif- 
ferent strategy is needed. 

In the present paper we apply the technique of 'eclipse 
mapping' (Horne 1985) for the study of the complex light 
curves of the active Algol W Cru. The binary system and 
new 7-color photometric observations in the Geneva sys- 
tem will be described in Sec. 2 and 3, respectively. In Sect. 
4 we describe the model and calculations, while Sect. 5 is 
devoted to the discussion of our results. 

2. Description of W Cru 

The binary star W Cru (HD 105998, CoD -58°4431) has 

the longest period among the group of W Serpentis bina- 
ries. Its exceptionally long period, P = 198 d, makes its 
a difficult object to observe. Only a single spectroscopic 
study is published (Woolf 1962), and long systematic pho- 
tometric observations are scarce. However, the situation 
was considerably improved after Plavec's (1984) appeal 
for new photometric observations. As a result two photo- 
electric data sets of W Cru were secured at the Auckland 
Observatory in 1984 and 1985 (Marino et al. 1988), and at 
the Nigel Observatory from 1985 to 1991 (Pazzi 1993). In 
1985 a systematic photoelectric campaign of W Cru was 
initiated by Z. Kviz at the Swiss Telescope on La Silla. As 
an initial results, using Geneva photometry, a new deter- 
mination of the primary minimum was communicated by 
Kviz & Rufener (1988). This aided proper phasing of the 
satellite observations. 

Woolf (1962) has observed W Cru spectroscopically 
throughout the orbital cycle and detected only the pres- 
ence of a less massive component with supergiant char- 
acteristics of the spectral type Gllab. Its RV curve is 
well defined and enable Woolf to derive the system's mass 
fimction /(m) = 5.83 Mq. Superimposed on the super- 
giant spectrum are hydrogen Balmer emission lines and 
displaced absorption lines; very typical for the stars of this 
class (c.f. Struve 1944). Woolf recognized that the deeper 
minimum in the light curve is due to an eclipse of the su- 
pergiant by an opaque body. Apparently, a more massive 
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star in this system is 'invisible' - Woolf drew similarity of 
W Cm to its more famous counterpart (3 Lyrae. 

While known for the other members of the group of 
these long-period interacting binaries, period change is not 
yet found for W Cru (c.f. Zola 1996). Asymmetry in the 
primary minimum, and somewhat erratic light curve from 
the cycle-to-cycle as is seen in our Fig. 3, are, at least, 
some of the reasons for unconclusive results regarding this 
issue. 

An important step in the understanding of this system 
was made by the independent studies of Zola (1996) and 
Daems (1998). Zola (1996) concentrated on the analysis of 
the light curves secured by Marino et al. (1988) and Pazzi 
(1993), respectively, while the study by Daems (1998) was 
more elaborated, and besides an independent modeling 
of unpublished Geneva 7-color photometry, also included 
some new spcictrophotonicitric and spectroscopic observa- 
tions. In particular, Daems was able to construct the spec- 
tral energy distribution of W Cru from lUE UV region up 
to IRAS FIR wavelengths. In both studies disk-model was 
employed, Zola used an cylindrical a-disk model, while 
Daems had applied an analytical torus-like disk. Further 
disciission of their results is left for Sect. 5 where they will 
be compared to the results of the present work. 

3. Photometric data and variability of the 
components 

Geneva 7-colour photometric measurements (Golay, 1980; 

Rufenerfe Nicolet, 1988) of W Cru were obtained from 
December 1984 to July 1989, using the Swiss 70 cm tele- 
scope at La Silla Observatory (ESO, Chile), equipped 
with the photoelectric photometer P7 (Burnet & Rufener, 
1979). During this period, 378 measurements of weight 
q > I have been obtained (see Rufener, 1988, for the 
definition of the weight q). The photometric reduction 
procedure was described by Rufener (1964, 1985); the 
photometric data in the Geneva system are collected in 
the General Catalogue (Rufener, 1988) and its up-to-date 
database (Burki, 2006). In the case of W Cru, the all-sky 
reduction was improved by the measurement of a neigh- 
bouring comparison star, HD 101021 (KlIII), a well mea- 
sured standard, with a confirmed stability: 440 measure- 
ments over a period of 20 years, mean V magnitude 5.1374 
with a standard deviation av = 0.0046. The same value 
of ay can be adopted for the global precision of the mea- 
sures on W Cru. The data on W Cru are listed in Table I 
available at the CDS. 

The magnitudes in each of the seven filters are ob- 
tained from the visual magnitude V and the six colour 
indices in the following manner : 

i = V-[V-B] + [i-B] (I) 

with i representing one of the seven filters U, B, V, Bi, 

B2,Vi, G. Remember that the Geneva [U-B] and [B-V] 
indices are not normalized to zero for an AOV star as it is 
the case for the Johnson UBV indices. 



In Fig. I light curves in all seven photometric pass 
bands of the Geneva system are shown. The following 
ephemeris is used for the calculations of the phases: Min 1 
(HJD)= 2440731.84 -I- 198.537 x E after Kohoutek (1988) 
and the elements are derived from all available minima 
found in the literature. The primary minimum, to which 
phase = 0.50 is assigned throughout in this paper, and 
which is caused by the eclipse of the visible component, be- 
came deeper for the shorter wavelengths. However, light 
curve in the U-hand has much larger intrinsic scatter - 
large changes in the width and depth of the primary min- 
imum are evident for the light curve in the C/-band than 
for the other pass bands. Unfortunately, phase coverage of 
the secondary minimum is much poorer than primary and 
cycle-to-cycle variations in the secondary minimum are 
not evident as for the primary minimum. In fact, proper 
phase coverage in the secondary minimum has been se- 
cured in a single cycle. The width of the secondary min- 
imum when an opaque disk is eclipsing its stellar com- 
panion is larger than the width of the primary minimum. 
At shorter wavelengths the secondary minimum becomes 
not only deeper but evidently narrower. Moreover, some 
other peculiarities in the light curves of W Cru are clearly 
seen: maximum following deeper minimum is higher than 
preceding, there are asymmetries in descending and as- 
cending branches of both, primary and secondary, minima. 
Bumps and/or humps in the light curves at the phases 
= 0.2, 0.3, 0.6 are very pronounced, and, moreover, vari- 
able from cycle-to-cycle. 

The color variations (bottom pannels in Fig. I) are in- 
dicative of temperature changes. In quadratures when the 
supergiant is seen from the side [B — V] is the largest. The 
amplitude of the changes in [B-V] index is about 0.25 mag. 
As is evident from Fig. I changes in [U — B] are twice as 
large then in [i? — F], i.e. about 0.50 mag. The variations 
are more difficult to interpret in [U — B]. Large and er- 
ratic changes in [U — B] around the primary minimum are 
present as consequence of a large activity in this system on 
time scale comparable to system's orbital period. It seems 
that the source of these rapid fluctuations, which has dif- 
ferent amplitudes from eycle-to-cycle is screened in the 
phase of the secondary minimum. But due to lack of more 
intensive phase coverage our conclusion is rather tentative 
at this point (see also Fig. 3). 

4. Eclipse-mapping: model and calculation 

4.1. An outline of the model 

The eclipse-mapping technique introduced by Horne 
(1985) in combination with Doppler tomography (Marsh 
& Horne 1988) became an essential tool for the study of 
the structure and mass-flow processes in cataclysmic vari- 
ables (CVs). In particular, eclipse-mapping is used to re- 
construct the surface brightness distribution of the accre- 
tion disks in CVs, and is often assumed to be an indirect 
imaging method. Initially, only the information contained 
in the shape of the eclipse was used in the assembly of 
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the map of the accretion disk surface brightness distri- 
bution. However, other sources of the radiation have been 
completely ignored, e.g. the secondary (often cool) compo- 
nent. Rutten (1998) has improved the technique with the 
development of 3D eclipse-mapping which fits the whole 
light-curve instead of just the eclipsed part. 

In the present study we have followed Rutten's (1998) 
approach. But, unlike the case of CVs in Algols here accre- 
tion disk is not prevailing source of radiation. Moreover, 
the hot component has a relatively large radius and its 
screening of the accretion disk should be taken into ac- 
count. 

Previous studies of the light curves of W Serpentis bi- 
naries have shown that hot components might be exposed 
to various degrees. Therefore, in the construction of the bi- 
nary model we have started from the general case in which 
all three main radiative sources are presented: hot mass- 
accreting component surrounded by an optically thick disk 
and cool mass-loosing component which (canonically) is 
filling its Roche lobe. 

The stellar surfaces arc approximated by equipoten- 
tial surfaces. These surfaces are determined by two pa- 
rameters: the mass ratio q and the filling factor S (for 
description of the gravitational potential of two rotating 
stars in binary system see e.g. Hilditch 2001). We approx- 
imated an equipotential surface with a grid of triangles. 
This was done to avoid ray tracing when calculating visi- 
bility. Geometry of an accretion disk in our model is given 
by an a-disk model, i.e. the disk is determined by the disk 
radius and an opening angle (}. 

The temperature distribution on stellar sTirfaces is 
parametrized by the effective temperatures of the compo- 
nents. This constant temperature distribution is modified 
by two effects: the gravity darkening and the reflection ef- 
fect. For the stellar atmosphere in radiative equilibrium 
the gravity darkening exponent is 0.25 (von Zeipel 1924) 
while for the convective atmosphere is set to the value 
0.08 according to Lucy (1968). Since we are using triangles 
rather than a grid of points it is assumed that the acceler- 
ation due to gravity on the stellar surface can be consid- 
ered constant along the triangle's surface for a sufficiently 
small triangle, and is equal to the acceleration in one of 
its vertices. Temperature distribution on a disk surface is 
computed in a simple power-law fashion. This distribu- 
tion is independent of the vertical coordinate, i.e. there 
is only radial distribution of the disk temperature. For 
a-disk model the temperature of disk is TooR~^/^. 

Although the geometry of an accretion disk is axially 
symmetric, its temperature distribution, in general, is not. 
In the first phase of the minimization (see Sec. 4.3) it is 
axially symmetric and is a power-law in radius. However, 
in the second phase of the minimization the requirement 
of the axisymmetry and radial dependence is relaxed and 
the temperature distribution is allowed to be essentially 
asymmetric. In all cases, the distribution is independent of 
the vertical coordinate, i.e. there is only radial and angular 
(with respect to the disk axis) distribution of the disk 



temperature. For a-disk model the temperature of disk is 
Tooi?-3/4. 

Radiation from one component influences the temper- 
ature distribution on the other component. To increase 
the speed of the calculations the reflection effect is com- 
puted in an approximate way. Stars are assumed to be 
spheres of constant effective temperature. For example, 
when calculating the effect of reflection of radiation from 
the secondary (component 2) on the primary (component 
1), a triangle whose one vertex has the position P(.x, y. z) 
is considered. The radiation dilution factor is defined l)y 
W = Ail/A-n: = AA2/A-kcP', where is the solid angle 
within which the inner hemisphere of the secondary com- 
ponent is seen. The hemisphere has an area A2, and d 
is the distance between the triangle and the center of the 
secondary component. A is the albedo of the primary com- 
ponent, A = Q means total absorption of radiation, A = 1 
means total reflection. 

4.2. Light curve synthesis 

The light curve is synthesized using a two-dimensional 
pixel grid onto which the image of binary system at a 
particular position is projected. This subsection describes 
how this is done in our code. 

One of the biggest difficulties in computing the total 
flux from the binary system at a given orbital phase and 
inclination is the problem of determining which parts of 
the system are visible and should be taken into ffux de- 
termination. One approach is to approximate each surface 
element by a single point, then trace the ray of light from 
that point to observer and check whether it intersects any 
other surface on its way. While this method is as accurate 
as the point-approximation, it is relatively slow since it 
involves three-dimensional lines. 

We have taken a different approach. The approxi- 
mation of stars and disk with a network of triangles is 
ideal when projecting a binary system image onto a two- 
dimensional pixel grid since a triangle can be filled using 
a simple algorithm (Bresenham 1965) which does not in- 
volve neither multiplication nor division. Triangles that 
are facing away from observer are not considered. The 
remaining triangles are sorted by distance and projected 
in descending order onto pixel grid. If the closer triangle 
partially covers the further triangle, it simply 'overwrites' 
the other triangle's pixels with its own pixels. What re- 
mains after all triangles have been projected is a two- 
dimensional image of the binary system. Thus, the visibil- 
ity problem is solved. 

Since only normalized light curves are used in our code, 
we define the fiux from a single triangle as the product of 
its surface and the blackbody intensity for the triangle 
temperature. At a given orbital phase the total fiux from 
all the visible pixels is divided by the total flux at phase 
0.25. 

We also include the possibility of having a source of 
this "third light," i.e. an unobscured component which 
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Fig. 2. In the upper panel the synthetic light curve (solid 
line) calculated for the optimized disk image and system's 
parameters compared to the observed light curve of W 
Cru (circles) in the V^-band of Geneva system. In the bot- 
tom panel is the same but for the ?7-band of the Geneva 
system. The light curves cover cycle 33 counting after the 
epoch in Kohoutek's ephemeris (see Sec. 3). 



provides a constant contribution to the light curve in all 
orbital phases. 

4.3. Minimization by Genetic algorithm 

Each synthetic model of the binary system with an accre- 
tion disk can be characterized by its light curve which is 
observed by an observer located at a large distance from 
the binary system's center. In order to deduce the parame- 
ters of the binary as well as the details of the accretion disk 
temperature distribution one needs to find that synthetic 
model whose light curve best matches the observed light 
curve. In practice this means using a minimization proce- 
dure on a set of parameters which describe the model. In 
this work we have employed the genetic algorithm (GA), 
see e.g. Charbonneau (1995), as the minimization proce- 
dure. GA is very efficient in case of a large number of 



parameters and especially in case of the image reconstruc- 
tion. 

The basic idea is to consider a set of N binary system 
models which are described by their basic stellar (mass 
ratio, Roche lobe filling factors, stellar effective tempera- 
tures, inclination) and accretion disk parameters (radius, 
height) as well as the accretion disk temperature distribu- 
tion (image). From the standpoint of a general GA, each 
model is simply a set of real numbers which we call organ- 
ism in the rest of this section. 

In general, we define the fitness Xi of organism as 
a continuous monotonous function of its parameters. In 
our particular case we define the fitness as the function of 
the light curve synthesized from the organisms parameters 
and the observed light curve. Since it is desirable that the 
fitness increases as the difference between the two light 
curves decreases, we define fitness as 



M 



{s^{j) - Sobs(j))' 



(2) 



where Sobs is the observed light curve containing M points 
and Si is the light curve synthesized from the organism i 
at exactly the same orbital phases present in Sobs- 

The GA starts with the random population of N or- 
ganisms which are sorted according to their fitness. The 
most-fit (better) half of the population is selected for 
reproduction, so that 7V/4 "parent" pairs are randomly 
selected (with better organisms being more probable to 
be selected than worse (less-fit) organisms) and then two 
"children" from each pair are created to replace the worse 
half of the population. This "reproduction" process is 
performed analogously to the similar process in biology, 
i.e. the organisms are divided into chromosomes and then 
the chromosomes of the children are obtained from the 
chromosomes of the parents through crossover process: 
first child gets one part of the chromosome from the first 
parent and the second part from the second parent, while 
the second child gets the first part from the second and 
the second part from the first parent. 

This process converges very quickly to the nearest lo- 
cal minimum. In order to avoid the "stagnation" of the 
population in a local minimum random "mutations" are 
introduced, such as "copy-errors" during crossover, ran- 
dom changes of the individual genes (numbers), as well as 
the global mutations such as the "asteroid hit": a large 
number of the population is radnomly re-initialized. 

Our fitting process consists of two steps: 

1. GA minimization of the binary system and accretion 
disk parameters without minimizing the accretion disk 
image (which is assumed to be a power-law function of 
radius) ; 

2. GA minimization of an accretion disk image while 
keeping the binary system's parameters constant. 

This splitting of our method is necessary due to the 
limited computational resources. A more general approach 
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Table 2. Absolute parameters for the components and 
an accretion disk in W Cru. Quantities are derived from 
solutions made from V band light curves. 



Parameter 


Zola 


This 




(1996) 


work 


Mass ratio, q 


0.16 


0.19 


Inclination, i 


87.8 


88.2 


Mass of gainer. Mi [Mq] 


7.82 


8.2 


Mass of loser, M2 [M©] 


1.25 


1.6 


Radius of loser, R2 [R© 


76 


71 


Disc outer radius, Rd [R©] 


132 


124 


Disc thickness, Zd [R©] 


14.5 


17 


Disc outer temp., Td,out [K] 


1500 


3600 



would involve an intrinsically three-dimensional disk ge- 
ometry and would then simultaneously minimize both the 
geometry and the temperature distribution. However, the 
number of parameters that describe disk geometry would 
be much larger than in our approach. Unfortunately, these 
parameters would only make a difference in a relatively 
small part of the light curve near the minima, i.e. the 
amount of deviation from the symmetric shape. Therefore, 
one would need to use a prohibitively large number of or- 
ganisms with a much larger set of parameters to minimize 
the light curve. Instead, we use our two-phase approach 
where the first phase can be thought of as the zeroth or- 
der phase (fitting of the light curve with axisymmetric 
disk temperature distribution and an optional hot spot), 
whereas the second phase can be thought of as the first or- 
der correction (asymmetries in the disk temperature dis- 
tribution introduce irregularities in the symmetric light 
curve minima). We find this approach more practical at 
this time, due to the abovementioned limited computa- 
tional resources. 

5. Binary parameters 

As shortly described in Sect. 1, both Zola (1996) and 
Daems (1998) have analyzed the light curves of W Cru 
with a disk model. Zola used mean light curves con- 
structed from Marino et al. (1988), and Pazzi (1993) 
photoelectric photometry, respectively. Daems study was 
based on then unpublished Geneva photometry. Both 
studies have arrived to a quite consistent set of the param- 
eters. While in our initial calculations we have started with 
the parameters which have cover very broad range in the 
parameter space, it was immediately clear that solutions 
will be in the narrow range around the values specified by 
Zola, and Daems, respectively. 

Since we are looking for the optimal set of the param- 
eters in a multi-dimensional parameter space it is recom- 
mended to fix as many parameters as possible. Therefore, 
we assigned theoretical values to gravity brightening, limb 
darkening, and albedo coefficients. The effective temper- 
ature of the G supergiant, only visible stellar component, 
has been derived by Daems (1998) from the spectral en- 
ergy distribution in the broad wavelength span, from lUE 



UV to optical and IR photometry obtained at ESO. He 
has arrived to T2^eS = 5500 K. 

Since W Cru is a single-lined SB and only the mass- 
function is known, mass ratio has to be derived from the 
light curves. In our, as in previous studies, a conserva- 
tive assumption that mass-loosing component is filling its 
critical Roche lobe, is assumed. Therefore, its size is deter- 
mined by mass ratio. This fact makes determination of the 
mass-ratio more easy through ellipticity effect. Otherwise, 
mass ratio will be poorly determined. 

Whatever initial set of parameters used, the solution 
has converged to the case with no primary visible at all. 
The, final search for the optimal set of the parameters was 
made with the remaining two binary parameters: the mass 
ratio q, inclination of the orbit i, and the following three 
disk parameters; disk outer temperature Td^out, angle Pd 
defining disk semithickness, and disk outer radius rd,out- 

The optimal solution, as defined by minimization 
with the Genetic Algorithm (see Sect. 4.3 for details) has 
given for the mass ratio g = 0.19, and inclination i = 88? 2. 
This should be compared to q = 0.16, and i = 87?8 de- 
rived by Zola, and q = (0.175; 0.160), and i = (88?2; 87?6) 
derived by Daems depending on cycle. 

In combination with the mass-function /(to) = 5.83 
R0 (Woolf 1962) for given mass ratio and inclination 
one can calculate the masses of the components. We get. 
Ml = 8.2 Mq, and M2 = 1.6 Mq. The sum of the masses, 
with the period which is known, gives the separation of 
the components A = 306 Rq. Since in the semi-detached 
configuration a relative size of the Roche-lobe filling com- 
ponent is uniquely determined by the system's mass ratio 
the radius of the supergiant in W Cru is R2 = 80 Rq . 

The optimal geometrical disk parameters are: = 
0.405, and 0d = 7?8 which translate to Rd = 124 Rq, 
and the semithickness at the outer disk edge Zd,out = 17 
Rq . The later parameter also defines an upper limit on the 
radius of the invisible component i?i < 17 Rq. This sup- 
ports the hypothesis that the primary component which 
is completely hidden by a thick accretion disk, is a mid-B 
MS star. Tabulation of Harmanec (1988) for the B-type 
MS star of mass about 8 Mq gives i? ~ 4 Rq, hidden 
well inside thick disk we have found, in particular for al- 
most edge-on view. Fig. 3 is showing time series plot of all 
available photoelectric data along synthesized light curve 
on the basis of our solution from Tab. 2. Intrinsic varia- 
tions, as discussed in Sect. 3, are clearly seen as due to 
variations in disk properties, both radiative and geomet- 
rical. 

We have attempted to compute solution for the U- 
band light curve. While keeping binary and stellar pa- 
rameters at the same values as in the F-band solution, 
only disk parameters were allowed to be adjusted. As ex- 
pected, disk parameters differ for different passbands with 
the disk more extended both in radial and vertical direc- 
tion and with a higher disk-rim temperature in the t/-band 
{Rd = 185 Rq, = 22 Rq, and Ta.edge = 4600 K). This 
is not surprising since both minima in ?7-band light curve 
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Fig. 3. Time series of V band photoelectric measurements of W Crucis, all available data sets are shown; filled circles 
stand for Geneva measurements (this work), open squares for Marino et al. (1988), open triangles for Pazzi (1993), 
and open circles for ESO (Manfroid et al. 1991; Sterken et al. 1993). Full line represent synthetic light curves from our 
solution which are derived from observations made in a single cycle no. 33 (aprox. HJD 2447100 - 2447400). Cycles 
are designated according the zero epoch in Kohoutek's (1988) ephemeris. 



are almost the same depth and slightly more broader than 
in y-band light curve (c.f. Fig. 2). 

Zola (1996) calculated a disk edge temperature using 
Smak's (1992) formula. Depending on the light curves used 
for the photometric solution he obtained Td,edge = 1500 
K, and 2000 K, respectively. Disk rim temperature derived 
in our modeling differs from Zola's estimates, and in par- 
ticular for the U-band is considerably larger. Even seen 
edge-on the disk is an important source of radiation; as is 
seen in Fig. 2 the minima are almost equal. It is an eclips- 
ing body producing a deep primary eclipse, and quite a 
long secondary eclipse without visible eclipse-shoulders. 



6. An image of disk 

An optimal set of the binary and disk parameters give 
only the first approximation to the light curve. Symmetric 
model of the binary can by no means give a fit to the light 
curves which exhibit a complex structure with asymmetric 
ascending and descending eclipse branches, unequal max- 
ima, and a number of bumps/humps in the light curve 
(Fig. 3). This was a serious problems faced by previous 
modelers as already described in Sect. 1. Once an optimal 
set of the parameters is found we can proceed with the 
second step in minimization (Sect. 4.3). In this step our 
minimization procedure has to recover the surface bright- 
ness distribution of an accretion disk. The disk is divided 
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Fig. 4. Reconstructed disk rim images shown in the phases 
0.0, 0.25, 0.30, 0.50, 0.75, and 0.80 (from top to bottom) 
in V photometric passband. Corresponding synthetic hght 
curves for the system with this accretion disk image is 
shown in upper pannel of the Fig. 2. Some features are 
identified which produces bumps in the hght curves at the 
phases 0.2, 0.3, and 0.65, respectively. 



in a grid of triangles. In looking for an optimal map of 
the disk brightness distribution the intensities of triangles 
could be varied. The optimal map is a reconstructed disk 
image. The information contained in the complete light 
curve is used, not only the parts covered in the eclipse (in 
the case of W Cru this will be the secondary eclipse) . 




Fig. 5. Reconstructed disk rim images shown in the phases 
of the minima in U photometric passband. Corresponding 
synthetic light curves for the system with this accretion 
disk image is shown in bottom pannel of the Fig. 2. Radial 
extension of the disk is larger for U passband than V pass- 
band. Some similarities of disk rim images in both pass- 
bands are seen for the primary (deeper) minimum, phase 
0.50 (upper pannel), but it is not so clear for the sec- 
ondary minimum, phase 0.0, which is considerably deeper 
inthe U passband light curve comparing to the depth in 
V passband (c.f. Fig. 2). 

In Fig. 4 reconstructed disk images are shown for sev- 
eral important phases in V spectral passband. Since in- 
clination is high we see almost a disk rim. The disk rim 
appears rather clumpy with nonuniform brightness distri- 
bution. It is difficult to speak on some structures but on 
the image in phase 0.0 when disk rim is visible in full ex- 
tension some kind of long 'filamentary' or 'stream' struc- 
tures might be recognized. There is an impression that 
they are slightly tilted to the orbital plane. Similar clumpy 
structures are revealed from [/-band light curve modeling 
(Fig. 5). Initial disk parameters are listed in Tab. 1, and 
an image reconstruction is obtained in the same procedure 
as was elaborated for the F-band. 

Hydrodynamical calculations of the gas flow in the 
close binaries have been performed focusing mostly on the 
study of accretion disks in cataclysmic variables. However, 
Bisikalo et al. (2000) have made calculations for the model 
which has resembled the long-period interacting binary /3 
Lyrae. So far, this is the only hydrodynamical simulations 
conducted for long-period interacting binaries. Some im- 
portant conclusions from their work are as follows: 1) the 
matter of the gas stream from the Roche-lobe filling com- 
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ponent is redistributed into two parts, one that forms a 
quasi-elliptic accretion disk, and the other which forms 
a circumstellar envelope around the mass-gaining compo- 
nent, 2) the gas stream approaches the disk tangentially 
and does not lead to the formation of the usually postu- 
lated 'hot spot', 3) in fact, the interaction between the gas 
from the circumstellar envelope and the stream results in 
the formation of an intensive shock wave, and faces the ac- 
creting star, and 4) a part of the gas from the circumstellar 
envelope is deflected away from the orbital plane after in- 
teraction with the incoming stream and leaves the system 
while the other one moves in the orbital plane, encircles 
the accretor and then collides with the stream. It is evi- 
dent that hydrodynamical simulations of mass flow in the 
interacting binaries like (3 Lyrae have shown a variety of 
the structures like accretion disk, circumstellar envelope, 
'hot line' and cool stream, halo. 

Very recent polarimetric study of W Serpentis, a pro- 
totype of this isolated class of the strongly interacting bi- 
nary stars, has been accomplished by Piirola et al. (2005). 
This system resembles many similarities with W Cru de- 
scribed in the present work (an asymmetric eclipse, large 
perturbations in the light curve also outside the eclipse, 
color variations in and outside the eclipse, etc.). Their de- 
tail modeling of the polarization measurements has shown 
complex structure of the circumstellar envelope in which 
the primary (hot) component is embedded. While photo- 
metric and spectroscopic studies have revealed the exis- 
tence of a geometrically and optically thick accretion disk 
around the hot component, no disk is seen in polarized 
light probably because multiple scattering and absorption 
effects in the optically thick medium reduce the polar- 
ization of the light. The authors have localized scatter- 
ing 'spot' which is probably associated with the optically 
thin polar regions. Similar effect in polarization behavior 
would be produced by a jet emerging perpendicular to the 
disk plane. It is obvious that Piirola et al. (2005) polari- 
metric study of W Ser emerged to the model similar to 
/3 Lyr as proposed by Harmanec et al. (1996) from an in- 
terferometric study, and supported by spectropolarimetric 
investigation by Hoffman et al. (1998). 

Our disk image reconstruction is based on a model that 
includes only an accretion disk and probably suffers from 
the absence of additional structures, in particular stream 
of the matter from cool companion in direction of more 
massive component and/or its, already, formed accretion 
disk. But pseudophotosphere of this optically thick disk 
which has completely hidden hot component of the binary 
system has undertaken its role in the emission of radiation. 

7. Conclusion 

Systematic photoelectric photometry of the long-period 
eclipsing binary W Cru in 7-color Geneva system has been 
carried out at the Swiss Telescope at La Silla from 1984- 
1989 and is, for the first time, presented in this paper. 
Several cycles of this almost 200 days long period binary 
system have been covered. Signatures of the still present 



activity in this system are obvious: maxima of the light 
curves are unequal, asymmetry is present in the eclipse 
branches, humps/bumps are present in the light curve and 
migrate from cycle-to-cycle, the width of the eclipses are 
varying. Previous research on this binary is corroborated: 
modeling of the light curve with an accretion disk model 
has revealed geometrically and optically thick accretion 
disk in which a more massive and presumably hotter com- 
ponent of mid-B type is completely hidden. 

We have used 3D eclipse-mapping technique to recon- 
struct an accretion disk image. Since the binary is viewed 
almost edge-on our map revealed an image of accretion 
disk rim only. Its nonuniform brightness distribution re- 
veals rather clumpy structure. We shortly discussed re- 
constructed disk image with hydrodynamical simulations 
of the matter flow in semi detached binaries which have 
shown that a number of features and structures might be 
formed in such processes. 

Acknowledgements. We acknowledge constructive comments 
from our referee Dr. Petr Harmanec which improved this paper. 
The observational part of this work has been partly supported 
by the Swiss National Science Foundation. KP acknowledges 
financial support from Croatian MZOS through research grant 
#0119254. 

References 

Andersen, J., Pavlovski, K., & Piirola, V. 1989, A&A, 215, 272 
Bisikalo, D. V., Harmanec, P., Boyarchuk, A. A., Kuznetsov, 

O. A., & Hadrava, P. 2000, A&A, 353, 1009 
Bresenham, J. E. 1965, IBM Systems Journal 4, 25 
Burki, G., 2006, Geneva photometric database, 

http: / / obswww.unige.ch / gcpd / phl3.html 
Burnet, M., & Rufener, M. 1979, A&A, 74, 54 
Charbonneau, P. 1995, ApJS, 101, 309 

Daems, K. 1998, PhD Thesis, Katholieke Universiteit Leuven 
Golay, M., 1980, Vistas in Astronomy 24, 141 
Harmanec, P. 1998, Bull. astr. Inst. Czechosl. 39, 329 
Harmanec, P. 2002, AN 323, 87 

Harmanec, P., Morand, F., Bonneau, D., et al. 1996, A&A, 
312, 879 

Hilditch, R. W. 2001, An Introduction to Close Binary Stars, 

Cambridge University Press, Cambridge 
Hoffman, J. L., Nordsiek, K. H., & Fox, G. K. 1998, AJ, 115, 

1576 

Home, K. 1985, MNRAS 213, 129 
Huang, S.-S. 1963, ApJ 138, 342 

Kviz, Z., & Rufener, F. 1988, Inf. Bull. Var. Stars No. 3158 
Knee, L. B. G., Scarfe, C. D., Mayor, M., Baldwin, B. W., & 

Meatheringham, S. J. 1986, A&A 168, 72 
Kohoutek, L. 1988, Inf. Bull. Var. Stars No. 3272 
LinneU, A. P. 2000, MNRAS 319, 255 
Lubow, S. H., & Shu, F. H. 1975, ApJ 198, 383 
Lucy, L. B. 1967 Z. Astrophys. 65, 89 

Manfroid, J., Sterken, C, Bruch, A., et al. 1991, A&AS, 87, 
481 

Marino, B. F., Walker, S. G., & Budding, E. 1988, Ap&SS 150, 
103 

Marsh, T. R., & Home, K. 1988, MNRAS 235, 269 

Mimica, P., & Pavlovski, K. 2003, in ASP Conf. Ser. 292, 405 



10 



K. Pavlovski, G. Burki & P. Mimica: Indirect Imaging of W Cru 



Pavlovski, K., & Kfiz, S. 1985, Bull. astr. Inst. Czechosl. 35, 

135 

Pavlovski, K., Burki, G., & Mimica, P. 2005, Ap&SS, 296, 417 
Pavlovski, K., Schneider, H., & Akan, M. C. 1992, A&A, 258, 
329 

Pazzi, L., 1993, Inf. Bull. Var. Stars No. 3847 

Piirola, V., Berdyugin, A., Mikkola, S., & Coyne, G. V., S.J. 

2005, ApJ, 632, 576 
Plavec, M. J. 1980, in lAU Symp. 88, 251 
Plavec, M. J. 1984, Inf. Bull. Var. Stars No. 2524 
Plavec, M. J., & Koch, R. H. 1978, Inf. Bull. Var. Stars No. 

1482 

Plavec, M. J., Weiland, J. L., & Koch, R. H. 1982, ApJ 256, 
206 

Rufener, F., 1964, Publ. Obs. Geneve, A, 664 
Rufener, F., 1985, iu lAU Symp. Ill, 253 
Rufener, F., 1988, Catalogue of Stars Measured in the Geneva 
Observatory Photometric System, 4e Edition, Obs. Geneve 
Rufener, F., Nicolet, B., 1988, A&A 206, 357 
Rutten, R. G. M. 1998, A&AS 127, 581 

Sterken, C, Manfroid, J., Anton, K., et al. 1993, A&AS, 102, 
79 

von Zeipel, H. 1924, MNRAS 84, 665 

Wilson, R. E., & Caldwell, D. L. 1978, ApJ 221, 917 

Woolf, N. J. 1962, MNRAS 123, 399 

Zola, S. 1991, Acta Astron. 41, 213 

Zola, S. 1996, A&A 308, 785 

Zola, S., & Ogloza, W. 2001, A&A 368, 932 

Zola, S., Hall, D. S., & Henry, G. W. 1994, A&A 285, 531 



